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SECTION  I 


INTRODUCTION 


If  the  gravity  field  of  a  homogeneous  spherical  earth  were  the  only 
force  acting  on  an  object,  the  classical  Keplerlan  orbital  ele¬ 
ments  describing  the  size,  shape  and  orientation  of  the  orbit  of 
the  object  would  remain  constant  In  time.  However,  the  presence 
of  various  perturbative  forces  causes  the  classical  (osculating) 
elements  to  vary  with  time.  These  variations  can  be 
divided  Into  three  categories:  secular,  short  period  and  long 
period.  The  short  period  variations  appear  to  be  functions  of 
the  position  of  the  object  In  Its  orbit.  The  long  period  varia¬ 
tions  appear  to  be  functions  of  the  position  of  perigee  In  space. 
Keplerlan  elements  from  which  the  short  and  long  period  varia¬ 
tions  have  been  removed  are  referred  to  herein  as  mean  elements. 

The  mean  elements  change  monotonically  with  time  and  the  rates 
of  these  changes  are  the  secular  variations. 

Analysts  making  decisions  concerning  orbits  of  various  objects 
must  be  aware  of  the  fact  that  Keplerlan  elements  vary  with  time. 
For  example.  If  comparisons  are  desired  among  orbit  determination 
solutions  at  various  epochs,  these  comparisons  must  be  made  In 
mean  elements.  If  classical  elements  are  used.  It  may  be  errone¬ 
ously  Inferred  that  a  maneuver  had  occurred  or  that  one  or  another 
solution  was  unreliable.  Since  mean  elements  are  well  behaved  In 
time,  mean  elements  at  any  epoch  can  be  computed  from  those  at  any 
other  epoch  if  the  secular  variations  are  known.  However,  it  must 
be  remembered  that  the  mean  elements  do  not  represent  the  actual 
position  of  an  object.  For  example.  If  observation  station  look 
angles  are  desired,  they  must  be  computed  from  osculating  ele¬ 
ments.  Hence,  It  is  necessary  to  be  able  to  convert  back  and 
forth  between  mean  and  osculating  elements. 

In  summary,  the  classical  elements  represent  the  actual  position 
and  velocity  of  the  object  but  are  poorly  behaved  In  time.  The 
mean  elements  do  not  represent  the  actual  position  and  velocity 
but  are  well  behaved  In  time.  Since  each  of  these  characteris¬ 
tics  Is  useful  In  orbit  support,  some  means  of  transformation 
from  one  to  the  other  Is  desired.  The  following  sections  give 
tne  equations  for  these  transformations  and  show  examples  of 
their  applications  to  both  theoretical  and  real  data. 
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The  Kepler ian  elements  used  in  this  discussion  are  defined  by: 

a  »  semi-major  axis 

e  ■  eccentricity 

1  ■  Inclination 

u)  «  argument  of  perigee 

n  «  right  ascension  of  ascending  node 

M  ■  mean  anomaly 

The  symbols  for  the  classical  elements  are  those  used  above,  while 
the  symbols  for  the  equivalent  mean  elements  are  those  above  In 
conjunction  with  a  bar  superscript  (e.g.  a). 

It  Is  not  the  purpose  of  this  report  to  derive  all  the  mathemati¬ 
cal  expressions  used  for  the  various  computations,  but  rather  to 
show  what  the  computational  procedures  are,  together  with  numeri¬ 
cal  results  obtained.  The  references  are  recommended  for  deriva¬ 
tions  of  the  techniques.  This  report  shows  examples  of  the  prac¬ 
tical  applications  of  these  techniques. 
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SECTION  II 


SECULAR  VARIATIONS 


When  using  mean  elements  for  orbit  support  it  is  necessary  to  have 
the  secular  variations  or  rates  of  change  of  these  elements  with 
time.  Often  when  mean  elements  are  obtained  from  an  outside  agency, 
the  secular  variations  are  included  in  the  transmitted  information. 
A  typical  example  of  this  case  is  the  Space  Defense  Center  (SDC) 
5-card  bulletin,  an  example  of  which  is  shown  in  Appendix  I.  This 
bulletin  has  been  recognized  as  a  useful  technique  for  transmitting 
orbital  information  and  now  enjoys  wide-spread  use  over  the  Air 
Force  Eastern  Test  Range  (AFETR) . 


When  secular  rates  are  not  furnished  with  the  mean  elements,  or 
when  the  mean  elements  have  been  derived  internally,  it  is  neces¬ 
sary  to  be  able  to  compute  these  rates  as  accurately  as  possible. 
This  can  be  accomplished  analytically  if  the  equations  are  known 
or  numerically  if  mean  elements  at  several  different  epochs  are 
available.  The  analytical  technique  will  be  discussed  here,  while 
the  numerical  is  discussed  in  Section  IV. 3.  The  analytical  secu¬ 
lar  rates  are  computed  as  functions  of  the  earth’s  second  zonal 
gravity  harmonic  coefficient,  J2,  and  atmospheric  drag. 


The  secular  variation  of  the  mean  mean  anomaly  is  the  mean  mean 
motion.  This  (perturbed  by  Kozai’s  factor)  is  given  by  (Ref.  1); 

"  ■  [  ^  [1  "  I  -  f  (II.l) 

where  p  is  the  mean  semi-latus  rectum  given  by: 


p  «  a  (l-e^) 


The  rates  of  change  of  the  mean  argument  of  perigee  and  right 
ascension  of  ascending  node  are  given  by  (Ref.  2  and  3): 

.  3n  J2  (5  cos^l  -1)  (II. 2) 

0)  ■  —  . — ■  ' 

452 

.  3n  J2  a2  cos  1 

n  - - - -  (11,3) 

2p2 

The  rate  of  change  of  the  mean  inclination  is  assumed  to  be  zero. 
All  element  rates  are  mean  element  rates. 


3 


(11.12) 


2 


1+e 


a 


where 


where 


ifer] 

**  16,667  revs/day 

A  -  0,  if  e  >  0.06 
A  «  4,  if  i  <  0.06  and  n  <,  16.204 
A  *  13,  if  e  <  0.06  and  n  >  16.204 


Thus,  the  mean  element  secular  variations  can  be  obtained  analytically 
for  whatever  use  may  be  desired  of  them.  The  capabilities  of  this 
technique  can  be  demonstrated  by  attempting  to  reproduce  the  element 
rates  apoearing  on  a  SDC  5-card  bulletin.  The  results ^of  such  a 
test  are  shown  in  Table  I.  The  six  moan  elements  a  d  a  were  input 
and  the  other  rates  were  computed  from  them.  Notice  the  difference 

between  SDC  and  analytical  values  for  ^  effects  should 

cause  negative  rather  than  positive  values. 


Analytical  secular  variations  are  used  for  trajectory  generation 
in  Section  IV. 2  where  it  is  shown  that  better  results  are  obtained 
^  H 

when  Y*  Y*  '~2*  ^  zero  for  that  particular  example.  It 

could  be  that  these  terms  are  necessary  for  use  with  objects  more 
significantly  affected  by  drag.  In  any  case,  it  is  probably  true 
that  the  empirical  expression  for  ^  could  be  improved.  Analytically 

2 

and  numerically  computed  secular  variations  are  compared  in 
Section  IV. 3.  Appendix  III  shows  the  FORTRAN  coding  for  computing 
all  the  secular  rates  described  above.  All  cases  studied  have  shown 
that  best  results  are  obtained  with  the  above-memtloned  acceleration 
terms  set  to  zero. 
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TABLE  I 


COMPARISONS  OF  SDC  AND  ANALYTICALLY  COMPUTED 
SECULAR  VARIATIONS 


Parameter 


SDC 

Object  //4483 


a  (earth  radii) 
e 

I  (deg) 
w  (deg) 

0  (deg) 

M  (deg) 
a  (ER/day) 

1  (ER/day2) 
e 

2 

• 

i  (deg/day) 
w  (deg/day) 

•j  (deg/day^) 

U  (deg/day) 
j  (deg/day^) 
n  (rev/day) 
j  (rev/day^) 

^  (rev/day^) 

•  •• 

|j(tev/day‘*) 


Bulletin  //28 

1.06351376 

.032704 

48.3932 

129.4386 

247.0671 

233.5949 

1.157451x10"3 

1.574609x10"^ 

1.0527x10"3 

2.8643x10"'^ 

0. 

4.84710 
8.8975xl0“3 
5.34388 
9. 8094x10" 3 
15.53805068 
1.2682872xl0"2 

0. 

0. 


Analytical 

Input 

Input 

Input 

Input 

Input 

Input 

Input 

-  3.198782x10"5 

-  1.0527x10“^ 

-  2. 9094x10“ 5 

0. 

4.84601 

8.8955x10“^ 

-  5.34265 

-  9. 807 2x10“ 3 
15.53797780 

1. 268280x1 0“2 

2.45174  xlO“'* 

0. 
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SECTION  III 


PERIODIC  VARIATIONS 


111.1  Introduction 

The  short  and  long  period  variations  In  the  classical 
Keplerlan  elements  can  be  approximated  as  functions  of  the 
earth's  second  and  third  zonal  gravity  harmonic  coeffi¬ 
cients.  Two  methods  for  computing  the  short  period 
variations  have  been  used  at  the  AFETR  Real  Time  Computer 
System  (RTCS),  one  developed  by  Y,  Kozal  (Ref  5)  and  the 
other  by  J.  B.  Frazer  (Ref  7).  The  Frazer  method  has 
also  been  used  to  compute  long  period  variations.  The 
two  methods  are  presented  below  along  with  discussions 
of  their  applicability. 

111. 2  Kozal  Periodic  Variations 

This  method  has  been  used  at  the  RTCS  to  compute  short 
period  variations  in  the  Keplerlan  elements  as  functions 
of  the  earth's  second  zonal  gravity  harmonic  and  Is  valid 
for  elliptical  orbits  only.  This  method  computes  the  short 
period  variations  from  the  mean  elements. 


First  the  eccentric  anomaly  is  computed  by  solving  Kepler's 
equation  Iteratively; 


E^ 

E 


M 

|1  _  -  eslnE^  -  M 

l-^cosE^ 


on  first  iteration 
on  succeeding  iterations 


‘  The  process  is  converged  when  two  successive  estimates  are 
within  some  epsilon  of  each  other.  Then  the  true  anomaly 
Is  given  by 

V  ■  2tan“^ 

The  mean  radius  from  the  center  of  the  earth  to  the  object 
Is  given  by! 

r  ■  a(l  -  e  cos  E) 

and  the  mean  seml-latus  rectum  Is  given  by: 
p  -  a(l  -  e^) 


PtecBiliwS 
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Kozal  applies  an  extra  perturbation  to  the  mean  semi-major 
axis  in  order  to  satisfy  the  following  relationship: 


[l  -  ^  [l  -  f  sln^l]  [l  -  if  ] 


where  a  ■  earth’s  semi-major  axis 
11^  -  earth’s  gravity  constant 

and  n  is  the  mean  mean  motion.  When  computing  osculating 
from  mean  elements,  first  the  extra  perturbation  is  removed 
then  the  short  period  variations  are  computed  and  added  to 
the  mean  elements.  The  perturbed  mean  semi-major  axis  is 
given  by: 


3J^a 


1  - 


--  -  f  Q  - 


(iii.i) 


Then  the  short  period  Keplerlan  element  variations  are 
given  by  (Ref  1): 


'“b  ■  T  {t  t  •  t  [(f)'-  '  ] 


/i\3  o-  -  - 

+  (^1  sin^i  cos  2  (v  +  w) 


(III. 2) 


4e 


f  Y  sin^I  cos  2  (v  +  al)^ 


] 


sin^I  cos  2  (v  +  0}) 

cos  2  (v  +  w)  +  e  cos  (v+2aj) 


+  j  e  cos  (3v  +  2(6) 


6i 


•) 

•^e  -  r 

r-j  sin2i  I  cos2  (vfw) 

p  L 

>] 


(III. 3) 


+  e  cos  (v  +  2a)  ) 


+  j  cos  (3v  +  2a)) 


(III. 4) 
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where  J  “  T 

e  lie. 

Notice  that  In  the  above  expressions  for  6a)g  and  6Mg  there 
are  eccentricity  divisors.  The  variations  in  these  elements 
are  undefined  for  circular  orbits  because  the  elements  them¬ 
selves  are  undefined. 

These  variations  are  added  to  the  mean  elements  to  yield 
quasi-osculating  (short  period  variations  only). 
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Appendix  IV  gives  the  FORTRAN  coding  for  this  algorithm 
and  Figures  1  through  6  show  theoretical  results  obtained 
by  the  method.  The  figures  show  the  classical  Kepler ian 
elements  and  quasi  mean  elements  (short  period  variations 
removed)  for  one  orbital  revolution  of  a  satellite.  An 
orbit  was  numerically  generated  with  initial  conditions 
in  the  classical  elements  given  in  Table  II,  with  no  drag 
and  using  the  gravity  model  shown  in  Table  III.  The 
Kozai  algorithm  was  applied  to  points  taken  from  the 
theoretical  trajectory  at  five  minute  intervals.  In 
this  case  osculating  elements  were  being  transformed  to 
quasi  mean  elements.  Since  the  variation  equations  are 
functions  of  the  mean  elements,  the  solution  was  per¬ 
formed  iteratively.  Secular  variations  in  the  elements 
were  removed  for  ease  of  plotting. 

Due  to  the  problems  of  this  algorithm  associated  with 
near  circular  orbits,  this  method  is  not  recommended  for 
general  orbit  support. 

The  effect  of  Kozai 's  extra  perturbation  can  be  seen  in 
Figure  1,  where  the  quasi  mean  is  displaced  from  the 
center  of  the  plot  of  the  osculating  semi-major  axis. 
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TABLE  II 

THEORETICAL  TRAJECTORY^ INITIAL  CONDITIONS 


Classical  Keplerlan  Elements 
a  «  22981800.  feet 


e  =  .02 

1  «  28.  :  deg 
0)  “  9.  deg 
0-115.  deg 
M  -  0.  ■  deg 


I 


i 
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TABLE  III 


THEORETICAL  TRAJECTORY  GRAVITY  MODEL 


a  = 
0 

6378.165  km 

398601.2  kmV 

sec^ 

ja 

m 

C  * 
nm 

S  * 

nm 

2 

0 

-1082.3 

0 

3 

0 

2.3 

0 

4 

0 

1.8 

0 

2 

1 

0 

0 

2 

2 

1.68 

-0.64 

3 

1 

1.77 

0.19 

3 

2 

0.29 

-0.03 

3 

3 

0.15 

0.14 

4 

1 

-0.57 

-0.46 

4 

2 

0.06 

0.26 

4 

3 

0.08 

-0.003 

4 

4 

-0.008 

0.006 

^Multiply  all  values  shown  by  10  ^ 
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Figure  1.  Kozal  Semi-Major  Axis  Short  Period  Variations 
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Figure  2.  Kozai  Eccentricity  Short  Period  Variations 
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Figure  3.  Kozai  Inclination  Short  Period  Variations 
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Figure  4,  Kozai  Argument  of  Perigee  Short  Period  Variations 
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T^igure  5.  Kozai  Right  Ascens  .1  of  Ascending  Node  Short 
Period  Variations 
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Figure  6.  Kozai  Mean  Anomaly  Short  Period  Variations 


III. 3  Frazer  Periodic  Variations 
III. 3.1  Introduction 

The  method  developed  by  J.  B,  Frazer  (Ref.  7)  does 
not  have  a  problem  with  circular  orbits  as  does 
the  Kozai  because  the  variations  are  computed  in 
Cartesian  elements.  This  method  is  valid  for 
circular  and  elliptical  orbits  only  and  as 
described  herein  is  invalid  for  equatorial  orbits 
due  to  sine  inclination  divisors. 

First,  the  mean  Keplerian  elements  are  transformed 
to  mean  Cartesian  elements  as  follows: 

r  =  r  U 
-► 

T  7 

r  "  r  U  +  rv  V 


The  osculating  Cartesian  system  is  defined  as  follows 
X,  y,  z  are  position  components  in  an  inertial  geo¬ 
centric  righthanded  system  with  x  and  y  in  the  plane 
of  the  mean  equator  of  date,  x  toward  the  mean  Vernal 
equinox  of  date,  and  z  toward  the  north  pole  ("mean” 
here  refers  to  the  absence  of  nutations). 


cos  u  cos  Q  -  sin  u  sin  Q  cos  I 
cos  u  sin  2  +  sin  u  cos  Q  cos  I 
sin  u  sin  i 


-sin  u  cos  Q  -  cos  u  sin  2  cos  i 
-sin  u  sin  +  cos  u  cos  cos  i 
cos  u  sin  1 


where  u  ■  mean  argument  of  latitude  «  w+v  and  where 
r  is  the  radius  from  the  earth’s  center  to  the  object 
given  by 


l+5cosV 

f  is  the  rate  of  change  of  r  given  by 

fu 


and  rv  is  the  product  of  the  radius  and  the  true 
anomaly  rate  of  change  given  by 

rv  - 

where  p  ■  mean  semi-latus  rectum  ■  a(l-e^) 


|j^(l+ecosv) 


The  perturbed  or  osculating  Cartesian  vector  due  to 
short  and  long  period  variations  is  given  by 

r  -  (r+6r)(U+6$) 

r  -  (H6f )  (U+6U)+(r^6rv)  (V+6V) 

or  the  variations  alone  are  given  by 

6r  -  6rWr6V 

6r  -  6tOf6rvV+?6$+t^6V 


Using 

W  =  - 

then 


W 

sinO  sinI 

X 

W 

B 

-cosO  sini 

y 

w 

cosi 

z 

6U  ■  (6u*f6ftcosi)  V+  (61  sinu  -6ftcosu  8inI)W 
6^  -  -(6u+6ftcosI)  (6i  cosu  +61^sinu  sini)W 
and  hence 

6r  -  6r^  +  r(6u+6ficosI)V 

+  r(61  sinu  -  6J2cosu  slni)W 


6r  «  [6f  -  rv(6u+6J2cosI)^  U 

+  [6r^  +  f(6u+6ncosI)3  V 

+  ^f(6i  sinu  -  6ncosu  sini) 

+  r^(6i  cosu  +  6S^sinu  sinl)j  W 

These  variations  are  composed  of  short  and  long 
period  components  as 

6r  -  6r  +  6rT 
s 

-  -V 

6f  -  6f  +  6f, 
s  L 

which  are  computed  separately  as  shown  in  the  follow¬ 
ing  sections.  Once  the  short  and  long  period  varia¬ 
tions  have  been  computed,  osculating  Cartesian  elements 
can  be  obtained.  The  osculating  Cartesian  elements 
can  then  be  transformed  into  osculating  Keplerian 
elements  by  standard  techniques. 
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III. 3. 


2  Frazer  Long  Period  Variations 

Although  Frazer  develops  the  variations  In  a  Cartesian 
coordinate  system.  It  Is  Interesting  to  note  the  forms 
of  the  equations  for  the  long  period  variations  In  the 
Keplerlan  elements  (Ref.  7),  For  all  zonal  harmonics, 
J„,  n>2; 


-  0 


6eL  ■  Sj^(l-e^)8lnl  cosX^ 
61^  *  “Sj^e  cosi  cosX^ 

j^5-2n-u 

vcos^I  I  slnXC 

sln^I  J  X 


izii 


lOcos^I 

l-5cos^i 


-  -  -  10  V 

=  -S  e  slnl  cosl  r - 7^  -  ' '  ~ 77 

L  n  l-5cos^T  sin'll 


slnXg 

X 


6Ml 

where 


3 

-S  aini  ^  slnX5 

n  e  X 


n-2 


S  .  ,  8(n-.l2;.  _ Ve _  ^ 

^  3*2  J2p^*^(l-5cos^l)  X,y,v 


I 


n-v  -  - 

— _vi-l  v-l_ 

(-1)  (n+v)I  e  sin  I 


^M+v  ^n^ 


\ vy  m  _ ^ _ _ _ _ 


A 

sln^^  a  -(-I)"?  gin  Xw,  X  even 
X>1 

^  cos  Xo),  X  odd 


2 

cosX^  ■  (-1)  cos  Xa),X  even 
«  (-1)  sin  XaSjX  odd 
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X  -  2,4, . . . , (n-2)  ,  n  even 
"  (n"2)  ,  n  odd 

U  -  X,X+2,...,  (n-2) 

V  *  X ,  X*^ 2 ,  •  •  • ,  n 

■  earth’s  semi-major  axis 

As  can  be  seen  in  the  above  equations,  there  is  no 
first  order  long  period  variation  in  the  semi-major 
axis.  Since  there  is  an  eccentricity  multiplier  in 
the  expressions  for  the  long  period  variations  in 
inclination  and  right  ascension  of  the  ascending 
node,  these  variations  are  insignificant  for  near¬ 
circular  orbits.  There  is  also  an  eccentricity 
multiplier  in  Sn  for  all  n  except  n-3(J3).  For 
near-circular  orbits,  then,  J3  is  the  only  zonal 
gravity  harmonic  coefficient  which  contributes  sig¬ 
nificantly  to  long  period  variations,  and  then  only 
in  eccentricity,  argument  of  perigee  and  mean  anomaly. 
But  argument  of  perigee,  mean  anomaly  and  their  varia¬ 
tions  (note  eccentricity  divisors)  are  undefined  for 
circular  orbits. 

In  the  above  equation  for  there  is  a  divisor  given 
as 

l-5cos^I 

This  expression  becomes  zero  when  inclination  is 
approximately  63.5  or  116.5  degrees.  This  is  known 
as  the  critical  inclination.  At  the  critical  inclina¬ 
tion  the  long  period  variations  in  the  orbital  ele¬ 
ments  are  undefined. 


After  considerable  manipulation,  Frazer  reduces  the 
long  period  variations  as  functions  of  the  zonal 
harmonic  coefficient  J3  to  the  following  in  the 
inertial  geocentric  Cartesian  elements  (Ref.  7): 


where 

“2 


6?^  ■  ra2  |slni  (l+ecosv)sinu  5 

+sinI(2+ecosv)cosu  ^  +  cosi  e  cosv  ^  (III, 10) 


slni(l+ecosv)cosu  U 


+  sinl(slnu  +  esinw)  V  +  cosi  e  sinv  w|  (III. 11) 


i  ^  Cl 

3  e 


2J2P 
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Notice  that  in  the  expressions  for  the  Cartesian 
variations  due  to  J3,  the  effect  of  the  critical 
inclination  does  not  appear. 


Appendix  V  gives  the  FORTRAN  coding  for  the  Frazer 
algorithm  (includes  short  period  variations)  and 
Figures  7  through  13  show  theoretical  results 
obtained  in  the  Keplerian  elements  by  applying  the 
long  period  variations  as  described  above.  The 
example  used  is  the  same  as  that  previously  given 
for  the  Kozai  method  in  Tables  II  and  III.  The 
mean  longitude  variations,  6Ll,  given  in  Figure  13 
are  defined  by: 

6L.  «*  6a)_  +  6M_ 

ij  li  li 

The  figures  show  mean  elements  and  quasi  mean  ele¬ 
ments  (short  period  variations  removed)  for  twenty 
days  of  satellite  motion.  Secular  variations  were 
removed  for  ease  of  plotting. 


Originally,  long  period  variations  due  to  J2  and 
were  included  in  the  computations,  but  numerical 
results  obtained  when  comparing  trajectory  genera¬ 
tions  via  mean  elements  versus  numerical  integration 
(Section  IV.  2)  were  much  poorer  than  those  obtained 
with  J2  and  effects  eliminated. 


If  the  computations  of  these  variations  are  restricted 
to  near-circular  orbits,  the  above  equations  (with 
e»0)  reduce  to; 


a.2  r  sini  (sinu  U  +  2  cosu  V) 


sini  (cosu  U  +  sinu  V) 


These  are  the  resulting  approximate  equations  even  when 
J2  and  are  included  in  the  full  equations. 

These  approximate  equations  yield  long  period  varia¬ 
tions  with  errors  varying  as  a  function  of  the  eccen¬ 
tricity.  Table  IV  shows  some  t3rpical  errors,  encountered 
at  various  eccentricities,  caused  by  using  the  approxi¬ 
mate  rather  than  the  full  equations.  Since  these  are 
just  isolated  examples,  they  must  not  by  any  means  be 
considered  as  maximum  errors  which  could  be  obtained. 

The  errors  given  are  the  differences  between  elements 
obtained  after  applying  the  equations  III. 10-11  and 
the  approximate  equations  above  to  an  input  vector. 
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Figure  7.  Frazer  Semi-Major  Axis  Long  Period  Variations 
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Figure  8.  Frazer  Eccentricity  Long  Period  Variations 
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Figure  9.  Frazer  Inclination  Long  Period  Variations 
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Frazer  Argument  of  Perigee  Long  Period  Variations 
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Figure  13.  Frazer  Mean  Longitude  Long  Period  Variations 
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TABLE  IV 


Examples  of  Errors  Caused  by  Using 
Approximate  Long  Period  Variations 


Eccentricity 

0.02 

0.125 

0.6 

(ft.) 

0. 

0. 

0. 

6y  (ft.) 

0. 

0. 

0. 

6*  (ft.) 

207. 

1926. 

7861. 

6x  (ft. /sec.) 

0.18 

1.59 

1.10 

6y  (ft. /sec.) 

0.04 

0.44 

0.19 

(ft. /sec.) 

0.32 

2.34 

3.84 

6a  (ft.) 

6e 

0.000004 

0.000111 

0.000010 

6l  (deg.) 

0.0002 

0.0023 

0.0096 

6u  (deg.) 

0.0233 

0.0421 

0.0275 

6ft  (deg.) 

0.0019 

0.0009 

0.0332 

6M  (deg.) 

0.0215 

0.0382 

0.0013 
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111.3.3  Short  Variations 


The  short  period  y'®o“^hI^nl/c^£fi- 

°cL2\rtS  1,  -u. 

is  (Ref.  7): 


ecosv 


(III. 12) 


6r  “  ap  |8in^icos2u  +  (l-3cos^I)|l-  3+;^cos7 

5t  .  -Oj^  |^2sin7i  (l+ecoBv)7  aln2u 

r  ilzii>+  .(1-»^co^3|(iii.13) 

+  (l-3co87i)  esinv  |-  2  l+d-e^)^ 


2cos2u  +  2ecos(2u-v) 


5rvg  -  |2cos2S  +  2icos(25-v) 

+ecosv  cos2u  1  (l+ecosv)  , 

J  .fa  -  -  2+(l-e^) 

-(l+ecosv)(l-3cos^I)  +  ecosv 
pZ-2  (esinv)  ^  'll 

^  2[i+(]-i^)’G 

61^  -  asini  cosi  |3[cos25  +  icos(25-v)J 
+  ecos(2u+v)^ 

651  •  -acosi  ^6(v-M  +  esinv) 

-3|^in2S  +  isin(25-v)j-isin(2S-(^)j 


(III. 14) 


(111.15) 


(III. 16) 
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6u 


-Scos^I)  (v-M) 

+4fl-6co8^i  +  ^  --^esinv 

L  l+(l-e2) J 

.(iWi) 

+2(5co8^i-2)  esin  (2u-v) 
+(7cos^I-l)8ln2u  +  2co8^I  e8in(2u-v)j 


(111.17) 


fhere 

a 


“i-“L5U 

the  difference  between  the  true  and  mean  anoma..les  needed 

to  conmute  6Q  and  6u  is  given  by: 

8  ® 

(v-M)  =  (v-1)  +  (E-M) 

where  E  l8  the  eccentric  anomaly  and  (Ref.  8): 


(v-E)  =  ain" 


r /esinv  \ 

/4-e^)Vl+  ecosv  u 

Ki+5cosv/ 

'  1+  (1+62)*^  ^ 

and 


,5  gx  = 

^e-m;  i+ecosv 


Chese  short  period  variations  (specifically 
and  6r^.)  Include  the  effect  of  Kozai  s  perturbed 
aemi-major  axis  (Ref.  1)  as  given  previously  in  Section 
III, 2. This  additional  variation  does  not  have  to  be 
accounted  for  separately. 
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Appendix  V  gives  the  FORTRAN  coding  for  the  Frazer 
algorithm  (Includes  long  period  variations)  and 
Figures  14  through  20  show  theoretical  results 
obtained  in  the  Keplerlan  elements  by  applying  the 
short  period  portion  as  described  above.  The 
example  used  is  the  same  as  that  previously  given 
for  the  Kozal  method  in  Tables  II  and  III.  The 
mean  longitude  variations,  6Lg,  given  in  Figure  20 
are  defined  by: 

6L«  -  6a)  +  6M 

8  8  S 

The  figures  show  quasi-mean  (short  period  varia¬ 
tions  removed)  and  classical  Keplerlan  elements 
for  one  orbital  revolution  of  a  satellite.  Secu¬ 
lar  variations  were  removed  for  ease  of  plotting. 

If  the  computations  of  these  variations  are  restricted 
to  near-circular  orbits,  the  above  equations  (with 
e*0)  reduce  to: 

6r  =  ap  sin^I  cos2u 
s  ^ 

6f  =  -2a-  sln^I  sln2u 
s  1 

6rVg  =  Q2sln^I  cos2u  -  (l-3cos^i)J 

6i  =  3asini  cosi  cos2u 
s 

60  =  3acosI  sln2u 

s 

6u  *  -  •?’  (7cos^I-l)sin2u 
s  z 

These  approximate  equations  yield  short  period 
variations  with  errors  varying  as  a  function  of 
the  eccentricity.  Table  V  shows  some  typical 
errors,  encountered  at  various  eccentricities, 
caused  by  using  the  approximate  rather  than  the 
full  equations.  Since  these  are  just  isolated 
examples,  they  must  not  by  any  means  be  considered 
as  maximum  errors  which  could  be  obtained.  The 
errors  given  are  the  differences  between  elements 
obtained  after  applying  the  equations  III. 12-17 
and  the  approximate  equations  above  to  an  input 
vector. 
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Figure  14.  Frazer  Semi-Major  Axis  Short  Period  Varlati 
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Figure  15.  Frazer  Eccentricity  Short  Period  Variations 
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Figure  18.  Frazer  Right  Ascension  of  Ascending  Node 
Short  Period  Variations 
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Figure  19,  Frazer  Mean  Anomaly  Short  Period  Variations 
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TABLE  V 

Examples  of  Errors  Caused  by  Using 
Approximate  Short  Period  Variations 


Eccentricity 


0.02 

0.125 

0.6 

6r  (ft.) 

98. 

385. 

3383. 

6f  (ft. /sec.) 

0.05 

1.29 

0.75 

6rv  (ft. /sec.) 

0.26 

0.42 

5.42 

61  (deg.) 

0.0002 

0.0002 

0.0006 

6ft  (deg.) 

0.0030 

0.0073 

0.0064 

6®  (deg.) 

0.0053 

0.0392 

0.0206 

«x  (ft.) 

628. 

3189. 

3612. 

6y  (ft.) 

633. 

577. 

5675. 

6z  (ft.) 

725. 

1090. 

788. 

4x  (ft. /sec.) 

0.04 

2.20 

3.46 

(ft. /sec.) 

1.12 

3.72 

10.33 

(ft. /sec.) 

0.62 

1.85 

2.97 

6a  (ft.)  ^ 

257. 

181. 

28717. 

6e 

0.000006 

0.000023 

0.0002611 

6u)  (deg.) 

0.0482 

0.0280 

0.0352  1 

6M  (deg.) 

0.0426 

0.0181 

0.0084 

III. 4  Inverse  Computations 

In  both  the  Kozal  and  Frazer  forms  of  the  variation  equations, 
the  variations  are  computed  as  functions  of  the  mean  elements. 
If  it  is  desired  to  transform  osculating  to  mean  elements  the 
equations  cannot  be  used  directly.  Instead,  they  are  used  in 
an  iterative  sense  as  described  below. 

First  an  approximate  set  of  mean  elements  is  obtained  (usually 
the  osculating  elements).  These  are  used  to  compute  variations 
and  approximate  osculating  elements  ^re  obtained  and  coiiq)ared 
to  the  known  osculating  elements.  The  differences  between  the 
two  sets  are  the  errors  in  the  computed  osculating  elements 
caused  by  errors  in  the  approximate  mean  elements.  These  mean 
elements  are  then  corrected  by  the  amounts  of  the  differences 
concluding  the  first  iteration.  This  process  is  repeated 
until  some  convergence  criteria  are  satisfied. 

When  using  the  Frazer  method  and  convergence  criteria  of  1.0 
feet  and  0,001  ft/sec,  convergence  is  almost  always  achieved 
in  three  iterations.  When  using  the  Kozai  method,  convergence 
is  sometimes  difficult  to  achieve  for  near-circular  orbits^ 
After  forty  iterations,  significant  errors  may  still  be  present 
in  argument  of  perigee  and  mean  anomaly.  However,  a  method 
involving  matrix  inversion  (Ref,  9)  can  be  used  to  yield  quick, 
solutions  for  the  Kozai  method. 

Appendix  VI  shows  the  FORTRAN  coding  which  performs  the  itera¬ 
tion  control  when  going  from  osculating  to  mean  elements  with 
the  Frazer  method.  Table  VI  shows  the  results  of  using  the  .. 
iterative  technique  to  compute  mean  from  osculating  elements 
using  the  Frazer  method.  Although  in  the  example  shown,  the 
long  period  variations  are  larger  than  the  short  period  in  the 
Cartesian  elements,  the  reverse  is  true  in  the  Keplerlan 
elements.  Also,  the  long  period  variations  in  the  Cartesian 
elements  make  no  significant  contribution  toward  variation 
in  the  semi-major  axis. 
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MEAN  KEPLERIAN  VECTOR  22952977.  0.018754  27.984143  18.085366  108.609094  299.836032 


SECTION  IV 


APPLICATIONS 


IV. 1  Int reduction 

The  preceding  sections  have  presented  the  mathematical  expres¬ 
sions  describing  variations  in  the  orbital  elements  and  have 
shown  numerical  examples  obtained  from  the  applications  of 
those  expressions.  The  purpose  of  this  section  is  to  show 
how  knowledge  of  these  variations  can  be  used  to  advantage 
in  orbital  support.  At  the  A¥£TR  RTCS.  two  major  uses  of 
these  variations  have  been  exploited.  These  two  uses  are 
described  below. 


Traiectory  Generation 


The  major  use  to  which  element  variations  have  been  applied 
at  the  RTCS  is  that  pertaining  to  trajectory  generation. 

This  is  based  upon  the  premise  that,  given  elements  at  epoch, 
elements  at  any  time  can  be  obtained  if  the  element  variations 
can  be  computed.  This  can  be  easily  accomplished  using  the 
mathematical  expressions  given  in  the  preceding  sections. 


For  example,  suppose  that  the  position  and  velocity  of  an 
object  is  given  in  the  osculating  geocentric  inertial  Car¬ 
tesian  system  at  epoch  and  the  position  and  velocity  is 
desired  in  the  same  system  at  some  other  time.  One  way 
to  accomplish  this  trajectory  generation  would  be  by  numeri¬ 
cal  integration  of  the  equations  of  motion  where  accelerations 
in  the  Cartesion  system  are  computed  based  on  models  of  the 
earth’s  gravity  and  atmosphere.  But  such  a  trajectory  gen¬ 
eration  can  be  accomplished  through  element  variations,  also. 
Such  a  generation  using  the  methods  presented  in  the  pre¬ 
ceding  sections  would  require  a  four-step  process.  First, 
mean  elements  would  be  obtained  from  the  osculating  (peri¬ 
odic  variations  removed)  using  inverse  computations  (Sec¬ 
tion  III. 4).  Second,  element  rates  (secular  variations) 
would  be  computed  (Section  II).  Third,  these  rates  would 
be  applied  to  the  mean  elements  at  epoch  to  obtain  the  mean 
elements  at  the  time  of  interest.  Fourth,  the  mean  elements 
at  the  time  of  interest  are  converted  to  osculating  (periodic 
variations  added)  thus  completing  the  trajectory  generation.. 
This  procedure  is  shown  pictorially  by; 


Preceding  page  blank 
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Mathematical  expressions  for  accomplishing  the  first,  second 
and  fourth  parts  above  have  been  presented.  The  third  part, 
comprising  the  update  of  the  mean  elements  from  epoch,  t  , 
to  time  of  interest,  t,  is  described  below.  Often,  the 
mean  elements  and  their  rates  at  epoch  are  given  in  the  form 
of  an  Air  Force  Space  Defense  Center  5-card  element  set 
(SDC  bulletin).  Information  contained  on  this  bulletin 
includes  the  following  mean  elements  and  their  time  deriva¬ 
tions  at  epoch  (see  Appendix  I): 

a  *  mean  semi-maior  axis 
o 

A  *  derivative  of  a 

Y  *  half  of  derivative  of  a 

e  »  mean  eccentricity 

A  =  derivative  of  e 

*2  =  half  of  v.itive  of  e 

I  =  mean  incxin?'.  Ion 
0 

i  =  derivative  of  I 

=  mean  argument  of  perigee 

u)  =  derivative  of  w 

Y  *  half  of  derivative  of  w 

=  mean  right  ascension  of  ascending  node 

=  derivative  of  Q 

Y  =  half  of  derivative  of 

=  mean  mean  anomaly 
=  mean  mean  motion  (derivative  of  M) 

=  half  of  derivative  of  n 


M 

o 

n 

n 

2 
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^  -  sixth  of  derivative  of  A 
o 

■  twenty  fourth  of  derivative  of  n 

As  can  be  seen,  these  bulletin  parameters  are  nothing  more 
than  the  coefficients  of  a  MacLaurin’s  series  expansion 
about  epoch.  The  mean  element  update  can  then  be  accomplished 
by; 

At  -  t-t 

o 

a  ■  a  +  aAt  +  ^At^ 

o  2 

e  *  e  +  eAt  +  “fAt^ 

o  2 

I  -  I  +  iAt 

o 

0)  ■  0)  +  wAt  +  -^t^ 

o  2 

U  -  +  fiAt  +  ^t2 

o  2 

M  »  M  +  nAt  +  ^At^  + 

O  /  D 

•  ^ 

In  practice  1  and  ^  are  usually  zero. 

The  above  method  for  updating  mean  eccentricity  can  sometimes 
cause  problems.  The  secular  variation  in  eccentricity  for 
most  objects  is  caused  principally  by  atmospheric  drag  which 
has  the  effect  of  circularizing  the  orbits.  In  other  words 
A  and  ^  ars  usually  both  negative.  When  e^  is  small,  e  can 
be  negative  for  some  values  of  At.  This  creates  problems 
when  trying  to  transform  these  mean  elements  into  osculating 
Cartesian  elements.  This  problem  can  be  circumvented  if  it 
is  assumed  that  the  eccentricity  decay  is  due  to  a  decrease 
in  apogee  height  and  that  perigee  height  remains  constant. 

This  is  usually  a  valid  assumption,  especially  if  the  ele¬ 
ment  rates  were  computed  based  on  equations  derived  from 
that  assumption  (Section  II).  Then; 

e  =  1  -  ^0  (IV. 7) 

a 

This  method  of  trajectory  generation  has  a  major  advantage 
over  numerical  Integration.  Since  it  is  not  a  step-by-step 
procedure,  the  vector  at  time  of  interest  can  be  immediately 
computed  with  no  regard  for  vectors  at  intermediate  times. 

When  At  is  large,  this  advantage  results  in  a  considerable 
savings  in  computation  time.  Being  a  step-by-step  procedure, 
numerical  integration  is  subject  to  round-off  and  truncation 
errors  which  can  be  fatal  when  At  is  large. 


(IV. 1) 
(IV. 2) 
(IV. 3) 
(IV. 4) 
(IV. 5) 
(IV. 6) 
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The  mean  element  update  method  of  trajectory  generation  has 
proved  to  be  sufficiently  accurate  for  satellite  acquisition 
at  AFETR  Miplr  radar  sites.  The  FORTRAN  coding  for  this  pro¬ 
cedure  is  given  In  Appendix  VII.  This  technique  was  compared 
against  numerical  Integration  for  two  cases:  with  and  without 
drag.  Table  VII  shows  the  Initial  conditions  used  for  both 
methods.  The  resulting  orbital  characteristics  are  also 
shown  as  well  as  the  ballistic  coefficient  used  with  the 
U.  S.  Standard  Atmosphere  1962  model  for  drag  computation. 

The  gravity  model  used  with  the  numerical  Integration  is 
given  In  Table  III.  The  central  gravity  term,  iip,  and  the 
zonal  coefficients,  J2  and  J3  were  also  used  for  the  mean 
element  computations. 

The  mean  elements  and  their  rates  at  epoch,  obtained  from  the 
initial  conditions  as  shown  in  Section  II,  are  given  in  Ta  jle 
VIII.  The  differences  between  the  trajectories  are  presented 
in  Intrack,  crosstrack  and  radial  components  in  Tables  IX  and 
X.  Differences  are  shown  at  ten  minute  intervals  for  the  first 
revolution,  the  last  revolution  of  the  first  day  and  the 
last  revolution  of  the  week.  The  differences  obtained 
In  these  examples  are  due  to  the  use  of  tesseral  harmonics 
in  the  Integrated  trajectory  and  to  inadequacies  Inherent  in 
the  mean  elements  trajectory  generation  scheme.  The 
differences  are  obtained  by  substracting  the  mean  element 
updated  trajectory  points  from  the  numerically  integrated 
trajectory.  Note  that  in  both  the  vacuum  and  drag  cases, 
the  radial  and  crosstrack  errors  oscillate  about  zero  with 
ever  increasing  amplitude  and  very  little  secular  variation, 
while  the  intrack  oscillates  about  an  ever-growing  mean  value. 

Drag  trajecjjory  generations  were  also  performed  with  zero 
values  for  and  The  resulting  trajectory  errors 

in  radial  and  intrack  components  were  considerably  smaller 
than  those  shown  in  Table  X.  For  example,  mean  differences  at 
the  end  of  one  week  were; 

radial  =  -2226  ft 

crosstrack  «  2179  ft 

intrack  «  -131289  ft 

For  some  reason,  these  terms  adversely  affect  the  accuracy 
of  this  drag  trajectory  generation.  It  could  be  that  they 
are  useful  for  objects  with  high. ballistic  coefficients  or 
in  extremely  low  orbits.  However,  all  cases  studied  have 
shown  that  optimum  results  are  obtained  with  the  above 
mentioned  parameters  set  to  zero. 
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TABLE  VII 


THEORETICAL  TRAJECTORY  INITIAL  CONDITIONS 


X  «  17837622.  ft. 
y  =  11170525.  ft. 

2  =  4559553.  ft. 

X  ^  —14197.  ft. /sec. 

y  =  17945.  ft. /sec. 

i  «  11635.  ft. /sec. 


Trajectory  Characteristics: 

h  100  nm 
P 

h  =  150  nm 
a 

p  =  89  rain. 


Ballistic  Coefficient: 

B  =  ^  =  0.01  ft2/lb 

w 
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TABLE  VIII 


THEORETICAL  MEAN  ELEMENTS  AND  RATES 

-  1.03625778  earth  radii  (ER) 

-  .006798 
»  29.999538  deg 
=  16.219243  deg 

-  9.998874  deg 

-  8.781831 

-  -2.843243xlO"‘‘  ER/day 

“  -4.020530x10“® 

-  -2.725106x10““ 

“  -3.853477x10“® 

-  12.086671  deg/day 

-  5.758724x10“  ^ 

<■  -7.612567  deg/day 
“  -0.003627  deg/day^ 

-  16.14925696  revs/day 
■  0.00332323  revs/day^ 

-  3.208824x10“® 


*These  rates  set  to  zero  for  vacuum  trajectory 
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table  IX 

VACUUM  TRAJECTORY  DIFFERENCES 


Prediction  Time 
-  (min) _ 

0 

10 

20 

30 

40 

50 

60 

70 

/  80 
(1  rev)  90 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

,  1430 

(1  day)  1440 


Radial 


Cross track 


Intrack 


0. 

16. 

-  3. 
25. 
78. 

108. 

99. 

-  4. 
“  63. 
-230. 


0. 

16. 

32. 

67. 

97. 

22. 

26. 

93. 

201. 

193. 


174. 

116. 

55. 

-  52. 
-122. 
-150, 

-  77. 

-  52. 

-  53. 

-  32. 


-  75. 
1829. 
2870. 
2470. 

926. 

-  890. 

-  2253. 

-  2612. 
-  1677. 

108. 


0. 

-  32. 
-104. 
-217. 
-299. 
-406. 
-512. 
-625. 
-542. 
-228. 

4418. 

4122. 

4012. 

4048. 

4185. 

4356. 

4532. 

4719. 

4807. 

4847. 


(1  week) 
Last  rev 


9990 

10000 

10010 

10020 

10030 

10040 

10050 

10060 

10070 

10080 

mean 


-290 

89. 

426. 

561. 

448. 

68. 

-322. 

-666. 

-765. 

-477. 

-  93. 


10113. 

18084. 

17421. 

8470. 

-  4440. 
-15140. 
-18589. 
-13065. 

-  1106. 
11436. 

1318. 


34607. 

34722. 

34394. 

33737. 

33057. 

32590. 

32601. 

33235. 

34174. 

35046. 

33816. 
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TABLE  X 


DRAG  TRAJECTORY  DIFFERENCES 
(POLYNOMIAL  UPDATE) 


Prediction  Time 
(min) _ 

0 

10 

20 

30 

40 

50 

60 

VO 

80 

(1  rev)  90 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

(1  day)  1440 

9990 

10000 

10010 

10020 

10030 

10040 

10050 

10060 

10070 

(1  week)  10080 

Last  rev  mean 


Prediction  Errors 


Radial 

Crosstra 

0. 

0. 

24. 

16. 

34. 

32. 

112. 

67. 

208. 

97. 

238. 

22. 

162. 

-  27. 

-  38. 

91. 

-  208. 

201. 

-  391. 

194. 

-  941. 

22. 

1126. 

1886. 

2972. 

2854. 

3668. 

2385. 

2904. 

817. 

968. 

-  967. 

-  1193. 

-  2266. 

-  2639. 

-  2556. 

-  2612. 

-  1574. 

-  1066. 

209. 

-72055. 

19431. 

-53853. 

14706. 

-41571. 

2075. 

-39677. 

-10114. 

-50072. 

-18215. 

-67623. 

-17533. 

-84648. 

-  8300. 

-92839. 

5071. 

-87841. 

16024. 

-72103. 

19268. 

-66228. 

2241. 

iftj _ 

Intrack 

0. 

3. 

41. 

172. 

367. 

664. 

951. 

-  nil. 

927. 

373. 

-  29273. 

-  29395. 

-  32253. 

-  36866. 

-  41387. 

-  43973. 

-  43544. 

-  40406. 

-  36347. 

-  33547. 

"1654791. 

-1659008. 

-1685803. 

-1723967. 

-1756987. 

-1769944. 

-1757641. 

-1727493. 

-1695574. 

-1678511. 

-1710972. 
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The  differences  shown  above  resulting  from  the  use  of  analyti¬ 
cal  secular  variations  are  considerably  larger  than  those 
obtained  from  numerical  computations  using  real  data  (Section 
IV. 3,  Table  XIV). 


There  is  good  agreement  between  numerical  and  mean  element 
trajectories  for  the  vacuum  case  (Table  IX).  Therefore  the 
large  intrack  differences  obtained  for  the  drag  case  must 
be  due  to  either  drag-affected  secular  rates  or  their  appli¬ 
cation  to  the  trajectory  generation  problem.  If  a  logarith¬ 
mic  rather  than  polynomial  update  scheme  is  used 
the  intrack  differences  are  considerably  decreased.  The 
logarithmic  scheme  is  based  upon  the  fact  that  the  mean 
mean  motion  is  a  logarithmic  function  of  time: 

• 

n  ■  n  +  log  (1+2  ^At) 
o  e  z 

and  the  mean  semi -major  axis  is  given  by: 


The  mean  mean  anomaly  can  then  be  updated  by: 


where 


and 


M  **  M  +  n  At 
o  a 


n  *  n  +  An 
a  o 


At 

An  =  log^ 


(1  +  y  t)  dt 


which,  after  changing  limits  and  applying  a  four-point  Gauss* 
method  of  approximate  quadratures,  becomes  (Ref.  10): 

A 

An  =  I  A^log^l  +  fAt  Xj 
1**1 


where  and  are  quadrature  values: 


1  0.3A785A85 

2  0.65214515 

3  0.65214515 

4  0.34785485 


0.93056815 

0.66999052 

0.33000947 

0.06943184 


Using  this  logarithmic  scheme  of  mean  mean  anomaly  update,  the 
differences  shown  in  Table  XI  were  obtained.  As  can  be  seen, 
the  Intrack  differences  are  much  improved  over  those  given  above 
which  suggests  that  this  method  should  be  extensively  studied 
both  experimentally  and  analytically.  However,  preliminary 
results  indicate  that  the  improvement  shown  is  typical  of  what 

can  be  achieved. 
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TABLE  XI 


DRAG  TRAJECTORY  DIFFERENCES 
(LOGARITHMIC  UPDATE) 


Prediction  Time  _ Prediction  Errors  (ft) 


(min) 

Radial 

Crosstrack 

Intrack 

0 

0. 

0. 

0. 

10 

24. 

16. 

3. 

20 

34. 

32. 

-  41. 

30 

112. 

67. 

-  172. 

40 

208. 

97. 

-  367. 

50 

238. 

22. 

-  664. 

60 

162. 

-  27. 

-  951. 

70 

-  38. 

91. 

-  nil. 

80 

-  208. 

201. 

-  926. 

(1  rev)  90 

-  391. 

194. 

-  373. 

1350 

-  959. 

22. 

-23497. 

1360 

1045. 

1885. 

-23532. 

1370 

2824. 

2852. 

-26214. 

1380 

3479. 

2384. 

-30588. 

1390 

2722. 

816. 

-34860. 

1400 

838. 

-  966. 

-37246. 

1410 

-  1250. 

-  2265: 

-36704. 

1420 

-  2638. 

-  2554. 

-33539. 

1430 

-  2596. 

-  1572. 

-29494. 

(1  day)  1440 

-  1088. 

209. 

-26677. 

9990 

-  6191. 

18850. 

61152. 

10000 

7573. 

14260. 

58903. 

10010 

16150. 

1992. 

39836. 

10020 

16693. 

-  9791. 

14265. 

10030 

7703. 

-17655. 

-  5041. 

10040 

-  6151. 

-17014. 

-  7541. 

10050 

-18414. 

-  8062. 

9340. 

10060 

-22944. 

4920. 

38732. 

10070 

-17505. 

15544. 

67476. 

(1  week)  10080 

-  4795. 

18673. 

82652. 

Last  rev  mean 

-  2788. 

2172. 

35977. 
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IV. 3  Orbit  Determination 


Another  use  to  which  element  variations  have  been  applied  at 
the  RTCS  is  that  pertaining  to  orbit  determination.  This 
technique  has  been  put  to  full  use  in  the  CDC  3100  com¬ 
puter  program,  CASS  (Ref  10) . 

Orbit  determination  using  mean  elements  is  accomplished  as 
follows.  First,  each  pass  of  tracking  data  is  reduced  to  a  . 
vector  (position  and  velocity)  using  some  method  of  orbit 
determination.  Second,  all  such  vectors  are  transformed  to 
mean  Keplerlan  elements  using  the  inverse  computations  des¬ 
cribed  in  Section  III. 4.  Third,  polynomials  in  time. are 
passed  in  a  least,  squares  sense  through  like  elements,  one 
from  each  pass.  The  polynomial  coefficients  thus  determined 
are  identical  in  definition  to  the  parameters  of  an-SDC  bulle¬ 
tin  (Appendix  I).  These  coefficients  then  represent  the  multi¬ 
pass  solution  and  can  be  used  to  predict  future  positions  of 
the  object  as  shown  in  Sec^on  IV. 2. 

In  CASS,  however,  the  polynomial  coeffj  tents  are  not  the  ones 
appearing  on  the  final  transmitted  bulletin.  It  is  necessary 
to  analytically  recompute  some  of  the  parameters  in  order  to 
preserve  known  relationships  among  them  (e.g.  a  and  4).  Some 
of  these  relationships  are  given  in  Section  II.  Examples  of 
actual  single-pass  mean  elements  used  and  the  polynomial  and 
analytical  bulletin^  thus  obtained  are  given  in  Tables  XII, 

XIII  and  XIV.  Prediction  errors  for  the  analytical  bulletin 
are  given  in  Table  XV  in  radial^  cross track  and  in track  com¬ 
ponents.  The  solution  and  prediction  errors  Were  obtained 
using. actual  tracking  data  from  AFETR  Mipir  radars.  These 
errors  are  caused  by  several  factors.  First,  the  mean  elements 
and  their  rates  at  epoch  contain  errors  due  to  their  computa¬ 
tions  being  based  on  imperfect  tracking  data.  Second,  the  vec¬ 
tors  against  which  the  predictions  are  compared  to  produce 
residuals  are  based  on  imperfect  tracking  data  and  contain ■ 
errors.  Third,  the  actual' forces  influencing  the  object  are 
not  completely .modeled  in  the  mean  elements  updating  scheme. 
Examples  of  some  not  included  in  the  computations  are  higher 
order  gravity  terms  and  atmospheric  density  variations.  Approxi 
mately  20.7  days  after  epoch  (17  August  1970)  a  large  geomag¬ 
netic  disturbance  occurred  which  significantly  effected  the 
prediction  errors  (Table  XV) . 

Notice  in  Table  XIII  that  the  single-pass  mean  element  resid- 
uies  from  the  polynomial  solution  show  periodic  errors  in 
mean  longitude  despite  (the  modeling  of  periodic  variations 
in  CASS.  The  cause  of  the  residual  long  period  variation  is 
not  known.  It  may  be  that  they  are  due  to  unmodeled  gravity 
or  drag  forces. 
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TABLE  XV 


CASS  ANALYTICAL  BULLETIN  PREDICTION  ERRORS 


Time  From  _ Prediction  Errors 


Radial 

Crosstrack 

3.2 

-  3. 

747. 

4.2 

466. 

-  2778. 

8.1 

-3206. 

3673. 

9.0 

-4083. 

3768. 

15.1 

-4344. 

2393. 

17.1 

-3757. 

775. 

21.2 

-1490. 

-13790. 
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-  8062. 

22.2 

-4287. 

-14080. 

24.1 

-3066. 

-  2328. 

25.0 
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2045. 

28.0 

-2452. 

-  9273. 

itrack 

4656. 

5892. 

8994. 

9482. 

8799. 

9588. 

25658. 

24480. 

31676. 

51068. 

67755. 

144416. 
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Comparisons  between  SDC  /f  “Xel. 

secular  variations  were  presented  in  Section  II,  Table  I. 

Here  note  the  differences  between  the 

ents  and  the  analytical  parameters  in  Table  XIV.„  Ttasg 
parameters  analytically  recomputed  were:  a,  7*  ®*  7* 

I,  and  n  using  methods  slmila-  to  those  given  In  Section  II 

(Ref  11). 
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(Key  on  following  page) 
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APPENDIX  II 


SUBROUTINE  DTPDA 


(Semi-Major  Axis  Decay  Rate  Computation  -  See  SECTION  II) 
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APPENDIX  III 


SUBEOUTINE  ELRAT 


(Mean  Elements  Secular  Rates  Computation  -  See  SECTION  II) 
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APPENDIX  IV 


SUBROUTINE  DELTA 

(Kozal  Short  Period  Element  Variations  -  See  SECTION  III. 2) 
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on  0000000000000. on. 


SuHRf)UTlN‘h  Dt-LTA  (  Afc ,  J'^ ,  f. ) 

THE^_  PURPOSE  OK  iHIS  SllRJ^Oi'T  T  NE  IS  TO  CO^Pt!  TF  SHPPT  PER.IOp 
ELEHENT  EVARiATlONS  USING  THE  KOZAI  MhTPOn. 

TnputVahe  '  '  '  ' 

APsEARTh  semi-major  A^IS  (FT)  . . 

J'->  =  EaPTh  gravity  SEf'UNI)  ZONAL  HARMOMC  COEFFICIENT 
M<l)=OH0li  mean  semi-major  AXIS  (I-T) 

Mrp)=MEAN  ECCENTRIC!  lY 

M(,\)=MEaN  I  NjCL_iNAT  1  nN_  (HEJS)  _  _ _ 

Mr4)=MEAN  argument  OF  PERIGEE  (OE'C) 

Mfs)=MEAN  RIGHT  aSceNSjon  OF  aSCENuING  NOPE  <nEC> 
Mf6i=MEAN  mean  anomaly  (DEG) 

Outputs  arl  differences  (UsculatIMg  -  mean>  Of  abovf  elements 
1*^  The  Same  Ugits  stqreu  in  D(.i)...I-HRU  U(6) _ _ 

DIMENSION  M(6),D(6) 
type  real  M,j2,JE 

DATA^PIsS  .l4lS926b.j6)  .  (F^=S/.^95779blJl  ) ,  ( OT=  .  333333333333  ) 
DATA  I  TTs  .66666^S666667  ) 

C  _  _ _ 

IPl  s  ^  PI 
JE=1 .5*J?*AE*AE 
E=H(?) 

IF 

1  E  =  n 

M(2)cE  _  _ 

GCf  1  o  '  ^  - 
?  I F  (  E  - 1 .  )  ^  /  3 »  A 

3  E=  .  999999y99\^99 
M(2)=E 

4  ESQsF^E 
OME2=l • -FSQ 

SOM='i";/0mF?*'^1 .5 .  . . .  . . .  ■" 

F  MsM<  <S)/FK 

9  IFIAHSF  (FN)-PI)  jO,lP,U 

c  Reduce  mean  aN"Maly  R/.nge  of  -pi  to  fPi 

13  FM  =  FM-.TPMaBSF(FM)/FM 
GO  rn  9 

compute  eccentric  'anomaly . .  ’  -  •  - . 

10  EPcFM 

DO  lA  J=l,10 

EAcER-IEP-MI^I^SINFCEPI-F  M)/(l.-M(2)*COSF  (FP)  ) 
IF(ARSF(EA-ER).LT.1,E"10^  go  to  20 
lA  EPsEa 

compute  true  ' A^ UMAI  Y 

20  V  =  ?.*aTAN2(S(;RTF  (  <1  .♦£  VU.-E)  )*SINF  (EA*  .5)  .COSF  (EA*  ,b)  ) 


IF(fM)  30  31.31 

30  FmsFm*TPI 

31  IF<V1  32.33.33 

32  VsV+TPl 

33_lF(Vj^TPl)  34.35.35  _  .  _  _ _ 

35  VsV-TPI 
34  ABsM<l) 

F leM(3)/FK 
FO«Mf 4)/FK 
FNoM(5)/FK 

COMPUTE  RaOIUS  and  SEMI-LATUS  kecTUM 

ir=A6*(  l  .  -M<T)iCO^Rt:A)  > 

P=AB*0ME2 

P2sP*P 

A0R3b( AB/R)**3 

compute  Various  FUFicTioNS  of  sin>-  and  cosinu  of 

C  inclination.  THUE  anomaly.  AI''D_ANGUMENT_  CF  perigee  _  _ _ 

'sf2s<SlNF<FI)‘)i*2 
SF 1=1 .-1.5*SI2 
SVsStNF(V) 

CVsCnSf (V) 

SO=SINF<FO) 

CO«CnSF(FO) 

■?v6='sv*co+cV«so‘' ■ 

CVOsf;7»CO-SV«SO 

S2V=P,*SV*CV 

C2V=2.*CV*CV-l. 

S3VsSV*<3.-4.*SV*SV) 

C3V=rv»(4 .*CV*CV-3,  ) 

STv^v*c3Tin::v;'53V .  . .  "  . . 

CAV=rv*C3V-SV*S3V 
S!5VsS2V*C3v*C2v*S3V 
C5V=r2V*C3V-S2V*S3V 
S20=p.*S0»CO 
020=2. *CO*CO-1. 

■■s'2Vo=2.‘»’svnVcvo . . 

C2VO»2.»CVO*i:vn-l. 

SV20«SV*C20*(:V»S20 
CV20=CV»C20-SV*S20 
S3V2nsS3V*C2l)*C3V*S20 
C3V2it  =  C_3V*C20_^S3V*S2n 

SvM2nsSv*r:20-cv*S20 ' 

S4V2nsS4V*C20*C4V*S20 
Si>V20  =  S5V»C20  +  C5V»S2n 
COMPUTE  kotai  short  periou  Variations 

[)(1)»( JE/AB)*(SFI*TT*( a0K3-S0M)*A0H3»3I2*C?V0) 

UEl»  IE*OME2/(  Ari»_Ad»E)»(OI  ♦SFI*(A0R3-S0M)  +  .‘i*A0R3*Sl2*C2V0) 
■  [JF:?  =  -ie»'5r2/'f2.»AB*E»P)*  <02''0*E*C'V2'0*0'T*EiC3V?0T 
0(2)=DEl-nE2 
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D(3>3jft»StNF(?.*r 1 )/(4.»fS)*<C2V0*H*CV20+0T*f*C3V20)*KK 
U01  =  '2.-V.5*Sl?)*(V-FM*t»SV>*SF  1*'  <l.-.2'5*(?SO)*SV/F« 

1  .^)*'!2V*h»S.SV/l2.  ) 

00?  =  -SV?n/t-«(  .25».St2*hStl*(  .5-.9375*312)  )♦.n6^5*e*SI2•SV^'^0 
f''0.3s-.,.‘>»(l.-2.'-*.Si2)*S?yu*S3(/2Q/e*(7.*SI2/12..-y.O*...  ...  . 

t  (1 .-19.*SI2/H. )/6,  ) 

O0'»=-3  75.S!2.S4V2O*.n6  •>‘j-fi*Sl2*S5V20-.375*SJ?*S20 
0(  4  )  a  JL-/P2*  <  n!n  +  L’U?*l)0  3  +  )*F'< 

0(3)s-JF:*C0SF  (FI  )7P2*(''-M*6*Sy-.5*S2tf0-.5«E*SV20- 
1  e*s.3v?n/6.  )»FK 

OMi  =  -SF  !•(  (j..  -  .^5*PSO)‘Sv  +  .5*6*S2V*(;SQ«S3V/i2.  ) 
dM2aM2*(  ;?5iSV2d*(l'.+i  .4<5*ESor-.o62b*FSO*SVH20“"  " 

1  7.«S3V20/12.» 

2  il.-eS-J/2fl.  )-.37‘>*E*SlV20*.n625*ESa*S5V?0+.375*F*S20> 
U(6)«  JE/<6*P2)*S0‘(TF(OM62)*(nMl  +  DM2)*FK 

tNn 


I 

1 


not  reproducible 
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APPENDIX  V 
SUBROUTINE  DELXYZ 


(Frazer  Long  and  Short  Period  Element 
Variations  -  See  SECTION  III. 3) 
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- "iiTTmijriNi;  -riHT'f YTT  *  t .  pi.j  .  DYTT  tS ; x7Z .  SSS')  '■  ■'  ■ 

i:  THJS  ACOH-Pr.i  MPaM  CARTFSIAN  •=l_rM(:’'jTS  ANI)_ 

‘"c  C0M3UTl='<'  ThP  n  if'FF  ifi'Ni:-:'?  -itlMFEM  'iFAi'j  A'Jn  OSc'Ji.AT  InG 

C  cartesian  SUliMtNTS.  '^^rgRtMCF  IS  FRaZHR.  SepT.  1966. 


c  iMPurs  aRu- _ ^ 

C  =  AXTs"  (FT)  '  '  . . 

n _ comstant  (Ff^^j/qFC**?) 

C  "  '  F"j23BAKTH  GRAVTrr  SFCWn  ZOMAu 'HA»?MniNlic 

C  FJ^s^AjUH  <5RAVITt  THHI)  /ONaL  HARMn^'nC  CO^rnCIFRT 

“  '  XY^(T/.  XYZ(?  )  i  Xy^'(i)  =  '1FA‘N  P'nSl  fTUM  r'OMpO-MFNTs  IN  All 

w  fcARTH-(:rE!\)T*i«<hU»i'JUAlORl  AL  ,  T  NFpT  I  ALi  R  I GHJ-H  ARDti) . 

C  '  TAfTn^  A V  ^  irirr G  A ,  ‘wlTTT  ('  1  r  "T n g  a*TP"  f H c  W I? N“A L 

C  _ _ egiJlNOX  ^NH  i6)  THRU_TmG  ^RTm  POCi  (FT) 

C"'  ■  *'  "  XY/(4yAT/r(‘5‘)>X'YZr6)=\/tL0CnY  'cOTiPn^^^^  Tn . 

C  A.U)VG  SYSTG:^  (FT/SEC) 

c"  ■  '■  . . . 


C  DjTPlJTS  ARE 

C  "  OPS  ( IT  I  Ou  I  :  'VVn  At  1  OnI?  ITvi  ■  CART1=<>T  AN  "EL^MEl^r  T 

C _  _  _^IJSj:^LArnMn-‘1EAR)  (FD  ANU  (Ft/SFC) 

C  ■  ■■  ■  . .  ’ 


x^^jS)j^0F'>v6) 

*  '  lF(FJ2.ffft.'T.  )  RJ  fo-io 

X=XY7( 1) 

Y  =  XYi:('^)  ■ 

_ 7=XYZ( 

xn=XY7(4)  . . .  ” 

_YP^XYZ(^J 

7n=xY/(6) .  . . . 

sHiJ=soRrr(F(-tu) 

compute  RaOIIIS . . 

RA0  =  30RTF  (X*X>VA-y  +  7*Z) 

COMPut?  IVv/FRTF  7fF''SPMi-MHJ0R  AXIS 

. .  a'I  N VT;  (  0  Vi  F MU-  R“ a  ;)■*  V ?  >  /  (  F  M  J ♦'«  A  0  ) 

|F(A(NV>LE,G,  )  \iO  Tj  iO 
HX  =  YVO-Yn^X  "  ‘  *  '  . . 


HYs-X*7i;)^^Xn*Z 
RZ  =  X*  YT-X')*  y" 
HXY23HX*TX^HY*HY 
H2sHXY'^~*-HZ*ri7 


not  reproducible 


_ 

C'OHPU  r  =  I  N'^L  fN  A  r  I  ON 

SINI_aSTRTF(HXY/!  J/R 

lF(STNi  .LT.<y.yo  317)  SO  TO  'l  l) 

COSICH7/R 

hsin'i  yVVs'iM  '  ■ 

cuMpurg  RT.  _ 

s'lNNsHX/HSlNl 


NOT  REPRODUCIBLE 
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COMMlj 


nuMPu 


COMPU 


COMPU 


nOMPli 


lOMPU 


COMPU 


rfrS:fa--4Y/Fr.TTTT - 

TF  or  LAjiTjnte 

S  I  Nija /‘/“(l^An  I  nT) 

"S  I  .  *S  ‘HU^COSU 

cos;?!i=i  .■;?.*srsij*si 
pi)or\/=x*xi]'*> 

15  PArr, 

RAT)l)5pnoT\/>RA’) 

_F  ^  i  i  V  =  R  IH  / r  M  0 

TF  StMi-LATUS  RRiJTii'l 

3aH^/FM(J _ 

SOP=SQRrt (P) 

pnos\/=(P/KAO)-i. 

■  OPf:CV  =  t  .•»'FC()SV"  " 

IF  TIM^_  fRUr 

R\/r)=SM'j+opPcv/sjp 
TF  ECCPNTPICITY  iOJA 
■■q^jspSpiV^FS  INV-^'tCOS 

FTAsSuorrCt .-P?i 
TF  TPDR  minus  MRAN  A 
S  V  M  P  s  (  F  sTn  V  /  0  Pi?  C  V  )  ■* 
ir< Aosr (SVMg) .nr.x, 
V  M  M  =  a“  sT^TfTsTNff )  rA 

TF  5H0RT_  PpRino  PL-RT 
“  RiT=<5i  ‘Mr*si'r;il 
r;i25Cosi*oosi 
sc'R^i 

PER=P/ aE 

"'ALf^HA  =  T;?5>*rj2/  c  P^E  * 
ALrAt_=ALrHAt»SMl)/SQ^" 

■  F S 2 u'^l '/ =  EC  1) S‘V  I  ‘Nl 2 U • 

FC2UM\/  =  pi;osVwCOS2U^ 
CS20P V  *s  I  M2U-* 
c(:2UPV3EcnsY^i:oS2U- 

1  <1.>ETA)  )  ) 

T  H  HP  a  I  j  P  EC  V  ^O'P'E'C  ^ 
_riE^<n=-  _  ALrAi*(  ? 

1'  (- JsVFTA^rElTP/Ti .+ 
nESRV03Sl2^^2.^C0SP 
nESRVi)af}irSRVn-S‘.TT2* 
1  .5»r(p?-2.*gS(NY*ES 
“  n  E  s  V  0  5  R  Y  n  *  A  L  P  H  A  *  n  p 
_nESl=Al-PHA-i3If\lI  *COS 
hEl>il5  “  aTFha“m!os  I  ♦  ( iS 
D  E  S  U  s  6  «  *  ( I ,  -  *  Cl  2  > 

OpStJanESU  +  SC  1  2* 

t  ES2'JMY 


<AD*HSInI  ) 


REPfj 


ODucn 


AiviOMALY  Rate 
XEII  _ 

tf*E::OSY'  . . . 

lil 

m')ma;.Y 

^ (eta*ooecv)/(i>fta))  ■ 

I  GO  TO  10 

. . . . 

jRaATIONS  IM  SppCiM.  FL:=*iPNTS 


=‘3IMY*C0S?U 

=  3rv!Y*S[.N2U  .  '  ' 

rSlNV^COSPO . 

=  31  *.iV*STiM2U 

^Ti ;  f  A7>TpE'fr:'+rc03v/ 


.*S!2*rFHP*9IN2U  +  SrI2^F3l  IJ* 

z'Ta)  )  )  '  ■  .  1- 

.i-^2 .  *cC?UhV-fEC0SV  ♦rnS^IJ) 

I  Ul+Tcos7Vl  (  ^\+Pf7)'/Ti7^rTA  ; 

l*MV)/(l.'*-PTA>  )  ’  V  ^ 

■>Rv'l 

1  3.  *(CU32il  +  i:C?0^./)-*'E.:2!JPV)  ‘ 

.  *(  V'1M-*-PSlNV;-:«5,*(^rt>l?iJ'^F32UMV)-ES2UP  /  ) 
*VM'1'^4.a(1.-6.*CI,?^,SCIP/(1,+ETA)  )'*4.SI  iM 
D  T  N  Y  *  F  c  n  3  Y7"(  i .  ♦  n  a  >  ?  V  *  ( •5".  r  12.-  2 . )  * . 
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- HFSua  I  ?-i  n**-!;  i  ’^^ij^rrvtrnri'ETprfFV 

nESj3  *,5*ALPHA* JFS  I  _ 

cof^Tnr'SHnKT  ^^RTj?<8AnaNS  im  Ca^f^Tan  «=l.SHfNT‘5i 

xu3cns‘i*CDS:^-siNU*SiNM*cnsi  _ _ 

YU=cn;s'!>^N^S!‘^«J»r:j^N*cosi  '  *  '  . . 

7UsSlN'J^Sl\l 

yV2-SlM*J*C0S'^-C0SUiDi\?*C0ST" 

YVS-^  ,MU*SII^N  +  C  )SU*  jOS>J*COSI 
7V  =  COSil*^rNT 

)(W  =  SIN!.!*SIM 

YH  =  -C0^'J*^TnI . .  ■ 

7W=Cns! 

T  E  H  P  5  n  «=  's  U  +  C i)  S  I  ♦  jFS  N 

TEM  =  S  I  MU’^DES  I  -*:.)SU* 3  I  M  I  •OESN 

nESX3l)PSR*XiJ+PA}*(TrHP*<Y+TE‘*1*XWr“ 

n£Sfsl)PSR*YU  +  RAJ*(TrMP*Y\/>TEM*Y»^) 

nYS7  =  nPSK*7J  +  PAjY(T;Yp*ZV>TEM*7'^)' 

TC  =  nESRl)-RVrj*rFMP 

TIJ  =  rjESRVl1-^R  Art-)*  *  "  "  -  --  -  . - 

TE  =  PA[)n*T(fM-^RV'l*(C03*J»nESI  i-S  lNU*SlNI*nPSN) 

nESYP  =  TC*YU  +  TO*  + 
n  E  TZT=  t  c  *Ti/>  T  T*  Z  V  "*•  T  r  7  4 

rOMPUTc  LO^G  PEHIOn  P.:RT»)XiUT  I ONS  IN  CARTESIAN  ELE'IE'^JTS 
'  AL?’AP=,5*Fj3/(rJ;?*Pr'ET  ■  ■*'  , 

_ cSO  =  Sp^J*ECOSV-COSlUESINV 

. .  . . . 

_ TH  =  MNJ_*  ULFA^^^ 

fC  =  Cn'ST*TALFA5?*:9C03  >  >■  '  . . 

_ TI)S“SIMI*0PECV*(  ALFa/^^^CJSJ)  _ _ _ 

TE=-.'il’MI  *(  alfa;?^(s1  xti  +  Esor)'  . . 

lF=_-r:oSJ*(  ALFA? ♦ESI  xV) 

“  nELX  =  RAiJ*(TAVx!j^T8*A'y^TC;*XN) 

_ nELrsRAD*(  rA*Yii^T8*  fV^rc*YR)  _  . .  _ 

nELZ  =  R  Af}*  (  T  A^Ztl  +  TR*/.  V  +  T3*ZW  ) 

_ T  EM  =  S:^i  \rj\^  P _ 

ni:LXn  =  TlfMw  (  rn^XJ+YF»  <V+ TFVx^jY .  .  "■ 

_r)EL Yn-TpH* ( r’2*Y j+TF*  Yv-f- rr^ry^) 

H  E  L  7  iV'  T  E  M  (‘t  13  ♦  z  I J  /  V  r  P  *  7  4  ) 

■30  rn  ?a 
1(3'  r^ONT  (  MUE 

[l6^^n(-‘SJ^-ni^Z=  3ESX  .i  =  DFSY(JrUE^7D=0  ._ 
nELX'Y3‘'LY  =  nF.LY=rjELXj’^OFUYn  =  (3EL7D  =  Ot  ‘  ’ 

_  ?0  f:9N T  IN-'E 

’'n¥s(  L7=!)irS<  +  l’3ELX  ' 

_ nES(?)=))ESY-*'jltLY 

"nES(X)3'i3FrsZ+^nELZ  ■ 

_ nES(4)5))(rSXl3>0FLXl3  _ 

nFS(5)=IJFSY;3-*-i)EirYL)'  . . 

_nFS(  6)rn£S_Zl3>Tr:LZl3 

. fnF .  ■  . . . 


not  reproducible 


APPENDIX  VI 
SUBKOUTINE  OSMCE 

(Iteration  Control  for  Computing  Mean  From 
Osculating  Elements  -  See  SECTION  III. 4) 


THjf  i .^►>os*p  ,)r'  ms  >iid-^jjriNh  IS  to  r^^«p«'JrF  •^eA^J 

OSi:Ul.A  f  ImJ'I  .:A-?n:ilA'l  *-.Lr^rNJS  ITFHaTI/cLY  o.ELxr/ 

.\i<^ 

A --f-A»<’r^'  SjfMi-MA  jV^  AXIS  (Ft)  . . . 

K’"j5tfAjn-l  COMSTaNT  (FT^^i/sCC**?) 

F  j;^=tA^TH '  I  r  f  sfconio  70nal  ha’^‘40^ic  cOErrlcicNr 
r  J,5  =  ifA»UH  <3-?WITf  HHO  ZONAL  CO'^FEiriENf 

xr/ (  i  )  .  <7z  ( >>)  *  Xir4(  .5)  =  OSCOLAT  INS  C.)MpONhNrS 

In  AN  r.ANTN-CfcNlFRPU,  FQUaTJN(Al*  IN‘=KTIAL^ 
NfSHT-HANDF  J’CA'NTESI  AN  SY.STE-1.*  wiTH  (1)  TO/<ArfO 
r.Hi^  VF'J^AL  rOUINOX  AND  (3)  HqU  TnF  NO’^TH  POL-  ( 
Xi^/(4)  .yYZ\S)  .XY/<6)  =  \/ElOCITY  C0rt®^^•E'ITS  IN 
A-IOVF  SYSTE1  iTT/SFO 


C _ O^iUTS  ANE  _  _ 

C~  '  <Y7*:P(T)-liS)V^EA'i  PJSlTinivj  AMi)  VEDriTY  roMPONENtS  IN 

C  3YSTEN  SIM!. AN  ID  aBOVF  (FT)  ^vD  (Fr/SFo) 

C  '  " 

n  I  I  ON  <  YZ  (  n  )  ,  |}P  >  (  6  )  ,  X  Y7M  (  6  >  .  OXYZ  (  A  ) 

C  ESTAH'.l'^H  IMTlAc.  tSilMATF  SF  MFaN  EL^  r-NTS 

yYZM(j)-XY.7rj)‘ 

1  OONflNNE 


1=0 

;?  co-viriviuE 
i  =  i«-i 

OAll  0‘^L<YZ(  AE.Fmu,  .  j;7,F  J3,  XYZM,  DES) 

ifTi  .\p.  Li'"sT'  ro  *4 

0  IMP.MVE  tSTIMATE  OF  -iEaN  ElFK-NTS 

no  ,5  J=l#6 

XYZM( j^=aY7( J)-0FS( J) 

.5  OONriNi’ifi 

_  so  ro 

4  c6Nn\"E 

iv^PNih/p  liSriMATr-  of  iEan  =^lfmfnts 

rO  •?  j=l»6 

nXfZ<J)  =  XY/U)-<YZMvJ)-tJ^=S(J) 

X  Y  Z  1  (  J )  =  X  Y  7  N  (  J  )  t )  X  Y  c  (  J  ) 

3  0O.Mj^J_N'it 

c . "Test  rno' c-j'wFNOtMCE 

iF<i,E''.iu)  no  ro 

no’T  j5i,i 

IM  AOSFOXYZCjn  .GT.l.  )  .',0  lO  o 
6  nn-NiriN'iE 
. J  =  4,6 

. l’F(Anbr(0XY7(J)  ).GT.'),n01)  SO  To'  2 


FD 
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APPENDIX  VII 


SUBROUTINE  UPDAT 


(Analytical  Trajectory  Generator  -  See  SECTION  IV. 2) 


Preceding  page  Wanli 
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«;il^RnijTTMtT  llPT),\T(T'^  VH,XN,  n  AY.  rn  ay,  A,P.  FI 
FmI-  PlfRf'O^b  OF  .rUS  I  mP  tS  TP  f  p^FORj!  thajl-ctuny 

PjCNFRaTuT'J  l:SlN(i  PEam  t.Lh'itiNT  SFCNlaW  vaRIa'^I'^NS  aNF'  a 

MACLa^'R  f  Nb  SERIFS  EM^/wSIO'Ni. 

ivpmrs  are 

■ . TTr‘’"r=i\rF.3"FR  pri^riov  OF  FPnCH  fDAvf;)’ 

r  i(''/=FRAcri()*ui  Rortiom  of  hRorn  rOAVs) 

X*ni)'-{A)  a:iJ  XMii)-(i;>)  aRp  The  '^paN  Ai^n 

iHPIb!  RATf:S  A!  b^OCH. 

. .  X'Ul)  =  *itANJ  SbMi-MAJilR  AXfC  (F'^') 

_^fUR)  =  ieA.'J  eccfwtricity 
<^nT)r1EAN  rNnUR*ATlO^'  VnPO'"  '  ‘ 

_XM(4)  =  'lbAN  aROOMEm]  OF  PpRiOE'-  (I)EG' 

’  X''t('5‘j,--»FANJ  RIGhT"~aSCE>fS!  o  !  OF 
_  •  ASCF.‘>»niruG  Nni.)F  OFr,) 

. ‘XM(61  =  -iEAM  -IhA'M  A^tUMAli^  rObG) 

_ JXJ'JI  Ijs^bAN  SbMi-MAJOR  AxtS  (E"?*CART!<  RADII/ 

'  XN('?)  =  AlToF=SbMV-MAjO  (FR>'DAY) 

X'^Jl  3)sa!)UO'.?  (F'</UAY**?} 

■  Xr.j(4j  =  ::0UT  =  FccbNfPlCITy  qAT^  (/OAY) 

XM(*5)r-n!)0'^  (/DAY**?) 

X  i(  6  )  -  jr)JT  =  ARG'jMEMT  (JF  PirP  1  OpP  R  aTP  («)fcO/nAY) 
_X'!(  7)- jOljO?  (ntO/nAY**^) 

X  g(  RVs ‘JODF- <inH"Y  ASCE\’5To  '  :>“F 
mOuE  UTb  (nE3/^»AY) 

'  X>i(R)'=Ml)002  (DrO/nAY**?) 

X  a  10  ) -•'i^'^  =  !*1rAM  MaTl»)fv|  (RpVS/D^Y^ 

X'J(ll)-'JUD2  (RbVS/iJAY*^?> 

12)  -i'iDD06  (REYS/DaY**-^) 

OA'YrlVfPJER*  PORT  1  ONI  O'  Ti^E  OF  ’JPDaTEO'  ^L^MEnTS  (DAYS) 
FOaYsE^aCT  (0  >IAL  -'IRTION  OF  T  I  M.P  JF 
ORnATFO  I-LE'iEMTS  (HAYS) 


O.lT'^l'TS  ARE 

A  =  tJ-MATEG  MPA.N  SrMI-M4J0R  AXIS  (t^?> 

bsU-’jVf!:;)"  r’paV  F^CtMTRIClTY 

P'l  -  IHOArEi)  MEAN  I- 'ICLINATION  (O^G) 

P.'ir  «Pi)A  TEn  'lEAN  aRv^'JmENT  OF  PERTo-p  (ObG) 

FjcIPOaFEO  "IbAN  '(IGHT  ASCENSION  OF  aSCENqInG  NODE  (ObO) 
rin«POATEil  '1E.AN  ibAN  AND'^ALY  (0^0) 

'nl'MEMStON'  TN(2)  ^XMC))  ,X\(12) 

nEUT  =  ( ^'AY  -  r  a  1 )  )> (F  ja Y  - T  a ? ) ) 
nEL  r  r  *nFLT 

A  =  X  vf  1  )  ^-DEi  T*X\’(2)  +  jFlT2*XM(:?) 

-I=XM(W 
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T7)-=XM‘(4'r^"(TT^DFLT^XN(7)*nEr:T?  . . . 

TN^XMth)  I-XN  (H)  *UbLT  ►A\(9)*nhLT? 

rM2XM  (  A  j  ♦.566  .’*  (  XN  ( 1  J  )  ^T4.xM(  1 3  )  *nbL  +  ( 12  )  ♦pF-LT  ♦OKI  F 2  ) 

RPDUc;^^  anomaly  iO  WA\GF  f)K  0  Ff:  .un 

PMSMIIOF  (r  M,36C  .  ) 

)  KM=FM+^^):i, 

*  PND 


not  reproducible 
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